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1 Introduction 

The Green functions G(x, x') of a partial differential operator, L : C°°(V(M)) -> C°°(V(M)), 
acting on fields <p G C°°(V(M)), i.e. smooth sections of a vector bundle V(M), over a 
(i-dimensional Riemannian manifold M with metric g is of great importance in mathemat- 
ical physics and quantum theory 0; H [§ EH- 111 particular, the singularities of the Green 
functions on the diagonal x = x' play a crucial role in the renormalization procedure of 
quantum field theory ||. 

In this paper we consider differential operators of higher orders of a special form, 
namely, the products of the second-order operators of Laplace type, i.e., 

L = F N ---F 2 F 1 , (1) 

where 

F t = -D +Q U i = 1, 2, . . . , N, (2) 

with □ = g M ^V M Vj, being the generalized Laplacian (or D'Alambertian, in hyperbolic 
case), V the covariant derivative on the vector bundle V(M), Qi = Qi(x) are arbitrary 
matrix- valued functions, i.e. endomorphisms of the vector bundle V(M), called the po- 
tential terms. Operators of this type arise in higher- derivative field theories (for example, 
in higher-derivative quantum gravity 0). We do not restrict the signature of the metric 
of the manifold, it can be Euclidean, signg = (+ • — h), leading to elliptic operators, as 
well as Minkowskian signg = ( — !-■■•+), corresponding to hyperbolic ones. 

The Green functions of the operator L are defined by requring them to satisfy the 
equation 

LG(x,x?) = 8(x,x'), (3) 

with L acting on the first argument of G(x, x') and 5(x, x') being the covariant delta- 
distribution on the vector bundle V(M). 
Here some remarks are to be made: 

i) First, in Euclidean (elliptic) case for noncomplete manifolds with boundary, the 
equation does not define a unique Green function and one has to impose some 
suitable boundary conditions. We assume that some appropriate boundary condi- 
tions are choosen, so that the Green function is defined uniquely. Since we are 
going to investigate only the singularities of G(x, x') near diagonal x = x', our anal- 
ysis will be purely local and the boundary conditions will not play any role in our 
consideration. 
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ii) Second, in pseudo-Euclidean (hyperbolic) case we will fix the Green function by 
the Wick rotation, i.e. the analytic continuation from the Euclidean case, which is 
equivalent to adding an infinitesimal negative imaginary part to all potential terms, 
Qj — > Qj — is. In other words, we consider the Feynman propagators. Thus in the 
following we deal mainly with the elliptic case — the formulas for the hyperbolic one 
are obtained just by analytic continuation. 

iii) Last, we will also assume for convenince the potential terms to be everywhere large 
enough, Qi(x) > const > 0, so that all operators Fi are elliptic and positive. Then 
the operator L is, of course, also positive. The singularities, do not depend, of 
course, on this assumption but it simplifies the intermediate computations. 

2 Heat kernel representation of the Green function 

The Green function of a product of operators (P is obviously expressed as the convolution 
of the Green functions Gj = F' 1 of each factor 

G = L- l = Ff 1 F 2 " 1 ■ • • = G 1 G 2 ■ ■ ■ G N . (4) 

Further, by using the heat kernel representation for the positive elliptic Laplace type 
operators @, |, |, [7J 

oo 

G,=F- l = fdtjUjitj), (5) 
o 

where 

UfiA^expi-tjFj), (6) 
one obtains the heat kernel representation for the Green function G 

oo oo 

G = Jdh-'-J dttfU^h) ■ ■ ■ U N (t N ), (7) 



or in the kernel form 

oo oo 

G(x,x') = jdtf-- J dtx J dvo\(yi) ■ ■ ■ J dvo\{yN~i)U\{ti\x } yi) 
o o M M 

xU 2 (t 2 \yi,y 2 ) ■ • ■U N {t N \y N _ l ,x'), (8) 



where dvo\(y) = dyg l l 2 (y) is the invariant volume element on M. 
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The behavior of the Green function G(x, x') near diagonal x — ► x' is determined by 
the behavior of the heat kernels Uj(tj\yj-i,yj) at tj — > 0, which is well known |3|, 0, ||. |7|. 
Thus, following J/J, we use a semi-classical ansatz for the heat kernel 

Wi\Vi-uVi) = {^r d/2 t- d ' 2 ^l\ yj - uyj ) exp ( n( % / '; //;) ) SlMVi-uVi)- ( 9 ) 

Here a(yj-i, yf) is the geodetic interval defined as one half the square of the geodesic 
distance between points yj-\ and yj and A(yj-i,yj) is the corresponding Van Fleck de- 
terminant [Q, [J. This ansatz reproduces the initial condition for the heat kernel, 

KmUjfolyj-!, yj) = %), (10) 

provided the function flj(tj\yj-i, yj), called the transfer function, has a regular limit as 
t - § g 

lim = ^(%-i> %)> (11) 

where V(yj-±,yj) is the parallel displacement operator along the geodesic. 
Now let us introduce new integration variables 

t = t 1 + --- + t N , &=j, (12) 
with the variables £j satisfying a constraint 

£i + --- + 6v = l. (13) 
The integration over the tj's takes the form 

oo oo oo ^ 1 1 

Jdh---J dt N f(tj) = Jdt — l - /#!•••/ d^ N 5(^ + ■ ■ ■ + ~ 1) /fa) (14) 

^ '' 

and, therefore, the Green function (^) can be written as 

G(*,*') = M- d/2 ^PzTjT Jdtt N ~ d ^e W (-%^) A(f|s,s') l (15) 

where 

A(t|a;,a; ) = (4tt) 1 ; 7 y dvol(yi) ■ ■ • y dvol^-i) A1/2 , ,x 

M M \ i ) 

oo oo 

x J d&... | ^(6 + "-- + ^-l)(ei---^)^ /2 



x exp (~2t $ ^' yi ' ■ ■ ■ n i(*Cik>yi) 

xn 2 (t&|3/i,2te) • ■■^ N (ti N \y N ^ u x'), (16) 
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where 

•<*,»,..., 6, 60 = ^ + + _ A (17) 

Remember, that locally, near diagonal, i.e. for x close to x', the Van Fleck determinant 
does not vanish A(x, x') ^ 0. 

3 Asymptotic expansions 



From the heat kernel representation (15) it is clear that the diagonal singularities of 
the Green function, when o~(x, x') — ► 0, are determined by the behavior of the function 
A(t\x,x') at small t. To study this behavior we remember, first of all, the well known 
asymptotic expansion of the transfer function near diagonal yj-i — > yj as t — > 0, || 

i ■ Vj ■• ~ E ( 18 ) 
fc>0 K - 

where a^\yj_i, yj) are the famous Hadamard-Minakshisundaram-De Witt-Seeley (HMDS)- 
coefficients ]1|, [II], 0, |11|]. They can be computed in form of covariant Taylor series and 
are known in general case up to |L2], ||, [13], |7|, |I4] . In particular cases, as in flat space, 
there results for higher-order coefficients. For a review of the methods for calculation of 
HMDS-coefficients and further references see [jlBl |16 |. 
By using the obvious equation 



: exp 



v{yji,yjY 



= S(y h , yj ) (19) 



it is not difficult to show that the function A(t\x, x') is regular in the limit t — > 

lim A{t\x, x') =V{x,x'). (20) 
Moreover, by changing the integration variables yj and using the asymptotic expansion of 



the heat kernel fll8|) one can obtain the asymptotic expansion of the function A(t\x,x') 



i.e. as t — > there holds, in particular, 

[(d+l)/2]-N-l ,_ xfc 

A(t\x, x') = Yl —n~c k (x, x') + R(t\x, x'), (21) 

where Ck{x,x') are some coefficients which are analytic in x — > x' and the rest term is of 
order 

R{t\x,x')^0(t^^- N ). (22) 



/. G. Avramidi: Singularities of products of Laplace type operators 



5 



The positivity of the operator L leads, in addition, to exponential falling-off of the function 
A(t\x, x') at the infinity, as t — > oo. Thus the integral over t in eq. ( |I5|) always converges 
at the upper limit t — > oo. 

Using eq. fl2ip one can calculate the integral (p~5|) over t and obtain an expansion of 
the Green function near diagonal, when o~(x,x') — > 0, 

G(*,*') = (An)-^ £ (-l) fc l * n , J A 1/2 (x,xOc fc (x,xO -^r 

+*(x,x'), (23) 

where 

¥(*,*') = M^ 2 ^^ ]dtt N - d ^e W (-^|^) i2(t|*,*'). (24) 

The first term in this formula is a polynomial in the inverse powers of o~(x, x') and 
gives the leading singularities of the Green function. Let us look at the second term, the 
function ty(x, x'). Off the diagonal, i.e. for x ^ x', or a ^ 0, the integral over t always 
converges at the limit t — > 0. However, on the diagonal, a = 0, the point is a bit more 
subtle. For odd dimensions, d, we have 

t N - d ' 2 - 1 R(t\x, x 1 ) ~ 0(r 1/2 ) (25) 



and, therefore, the integral ( pip over t converges at the limit £ — ► also on diagonal 
x = x', meaning that the function ty(x,x') is a regular smooth function on the diagonal, 
actually it is analytic in x close to x' and, hence, can be expanded in a covariant Taylor 
expansion. 

For even dimensions, d, in contrary, there holds 

t JV-«»/2-i i2 ( f | a . jX /)^ ( f -i) j (26) 

and, therefore, the integral (p4|) does not converge on the diagonal, for a = 0. It is easy 
to see, that in this case the function *S?(x, x') exhibits a logarithmic divergence 

t/ a (4vr)- d / 2 (-l) d / 2 - Ar + 1 Al/2/ . t , , , 

= (AT - 1)\T(* -N+l) ^ [-^^J+^re^x'^), 

(27) 

where ty ieg (x, x'\ fx) is a regular analytic function in x — > x' and \i is an arbitrary dimen- 
sionful parameter introduced to preserve dimensions. 
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Thus we obtained the complete structure of the diagonal singularities of the Green 
function of the product of N Laplace type operators 

\^]-N-l wrf Ar , n N i-N-k 



L 2 J -n(d _ at _ u\ 

G(x,x') = (An)- d ^ J2 (-1)' ° * } ^ /2 Mc k {x,x>) 



k=0 



d 
2_ 

fc!(JV-l)! ~ v ~'~ y ~ KV ~'~ ' \a(x,x') 



where 



^^ (N-t^-N+l) ^ ^ + Greg(X ' M) ' (28) 



' l)d/2 N+i Cd ( x x '\ for even d 
XM = { (29) 



0, for odd d 



Therefrom, we see that: 



i) the structure of singularities is determined by the coefficients Cf-{x, x') of the asymp- 
totic expansion of the function A(t\x,x'); 

ii) in odd dimenions there is no logarithmic singularity and the regular part of the 
Green function G reg does not depend on \x\ 

iii) for iV > d/2 there are no singularities at all and the Green function is regular 
analytic function on the diagonal, i.e. there exists finite coincidence limit G(x,x); 

iv) for N = [(d+ l)/2] there are no singularities for odd dimension d = 2N — 1 and 
there is only logarithmic singularity in even dimension d = 2N. 

The logarithmic singularity of the Green function is very important. On the one hand 
it determines, as usual, the renormalization properties of the regular part of the Green 
function. Indeed, noting that G does not depend on \x we obtain from (f2~8|) 

/i^G res = t r. — 72 rA 1/2 (x, x') X (x, x'). (30) 

^d\i g (iV-l)ir(f -iV + 1) v ' x 1 K ' 

On the other hand, the absence of the logarithmic singularity is a necessary condition 
for the validity of the Huygence principle for hyperbolic operators [17, 18, 19 1. 
The problem is now to calculate the coefficients C}~(x, x'). 



4 Particular cases 
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4.1 iV-th power of a Laplace type operator 

The simplest case is, of course, the case of equal potential terms, i.e. 

Q 1 =Q 2 = --- = Qn = Q- (31) 
In this case all operators Fj are equal to each other 

F a = • • • = F N = F = - □ +Q, (32) 
and the operator L is the iV-th power of the Laplace type operator 

L = F N . (33) 

Using the heat kernel representation we find 

1 °° 

G = F~ N = — — — ( dtt N ^U{t), (34) 



where U(t) = exp(-tF). Rewriting this in kernel form and comparing with eq. ( |15"D we 
find that the function A(t\x,x') in this case is just the transfer function determined by 
the ansatz @ for the heat kernel of the Laplace type operator F: 

A{t\x,x') = tt(t\x,x'). (35) 

This means, in particular, that the transfer function satisfies a nontrivial identity 

n /,| f\ fA s-(N-l)d/2 { A i/ x fr U ,A 1/2 (x, J/i) • • • A 1/2 (?/Af_i, x') 

Vt(t\x,x) = (4vr) 1 >' J dvol(yi) • • • J dvol^jv-i) AV 2 (x x') 

oo oo 

x J J d£ N 6(Z 1 + ■■■ + £„- l)^!---^)-*' 2 



x exp f-i$(ar, y 1 ,...,x'; £ u £ N )\ y%) 

xn(t6|2/i,2/2) ■ • •^(^7v|z/7v-i,a; / ) ) (36) 

where yi, . . . , x'; £i, £/v) is defined by eq. fll?]). From eq. ( |35| ) we find that the 
coefficients are equal to the usual heat kernel coefficients 

c k = a k . (37) 

Thus the extent to which the coefficients c k differ from the HMDS-coefficients a k ^ in 
general case is only due to the different potential terms. We see that the function A(t) is 
an 'averaged' (or weighted) transfer function of the Laplace type operators. 
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4.2 Covariantly constant perturbations of the potential terms 

Let us consider now the case when the potential terms Qj differ from each other by 
covariantly constant terms Xj which commute with all other matrices (the potential terms, 
the bundle curvature 1Z^ V = [V M , V„] and their derivatives), i.e., 

+ Xj, V^Xj = 0, (38) 

[\j(x),Q(v)] = M X )>KM] = [Hx),Xj(y)] = 0. (39) 
In this case we have 

Fi = F + X h (40) 

where 

F = -n+Q. (41) 

and, therefore, 

L= {F + X N )---(F + X 2 )(F + X 1 ). (42) 

The operators F{, even if different, commute and from the heat kernel representation (|7|) 
by using the change of variables flllj) we obtain 



oo 

G = L- l = J dtt N - l u{t) exp(-tF), (43) 

where 

^ g-t(Ai-Ajv) _ -y e -t(A 2 -A JV ) _ ^ g— *(^JV-l _ ^Jv) _ \ 

t(X\ — Xn) t(X 2 — Xn) £(Ajy_i — Xn) 

Rewriting this equation in the kernel form and comparing with eq. (p~5|) we find that the 
function A(t) differs from the transfer function Q(t\x,x') for the Laplace type operator F 
(p|) just by a prefactor u(t): 

A(t\x,x') =u(t)tt(t\x,x'). (45) 
The function u(t) is analytic in t. Expanding it in the Taylor series 

"(t) = E [ -J-"n (46) 

we find the coefficients 

\u} n a k _ n , (47) 



0<n<fc \'V 



where the HMDS-coefficients. 
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5 Algebraic approach 

As we have seen, the function x') is very closely connected with the transfer function 

Q(t\x, x') of an Laplace type operator. Moreover, in particular cases considered above the 
coefficients Ck are expressed linearly in terms of the usual HMDS-coefficients a^. Let us 
show now how one can compute the coefficients c& in general case. 
Let us define some auxilliary operators Zj, (j = 1, 2, . . . , N), by 



Zi 



F\G — F 2 1 • ■ ■ F N X 



Z 2 — F 2 Z 1 — F 2 F 1 G — F 3 1 • • • Fj^ 1 



Zn-i — Fn-iZn-2 — F N X 
Z N = F N Z N _i = 1 



(48) 



Now let us introduce a column-vector built from the auxiliary operators 

f g ^ 

Z 2 



(49) 



V 



The Green function we are looking for is the first component of this vector and can be 
obtained by multiplying with the row- vector U\ = (1, 0, . . . , 0) 



G = U\Z. 



Further, let us introduce a matrix F formed by the operators Fi 



(50) 



/ 



Fi -1 ... 
F 2 -1 ... 



\ 



(51) 



\ ... F N J 
Then the equations ( f48l) can be rewritten in a matrix form 



FZ = LLv, 



(52) 
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where 



n 



N 



(53) 





v 1 ; 



The operator F has the form 



F = -l-n+Q, 

where 1 is the unity N x N matrix and Q is defined by 



(54) 



Q 



( Qx -1 ... ^ 
Q 2 -1 ... 



(55) 



\ ... Q N J 

It is of Laplace type and is, therefore, nondegenerate, i.e. there exists a well defined Green 
function G = F~ l . Thus from the eq. (|5^) we have 

Z = GU N . (56) 

Finally, from (|50|) we get the Green function 

G = Tl\Z = U\GU N = tr (UG), (57) 

where II is a projector of the form 



n = n 



N 



o o 







... 

1 ... 



(58) 



and 'tr' denotes the usual matrix trace. Note that this trace has nothing to do with the 
trace over the bundle indices which are left intact. 

This shows that the problem is reduced to a Laplace type operator F with additional 
matrix structure. Using the usual heat kernel representtion for the Green function of the 
Laplace type operator we have from eq. ( |5TD 

a(x, x') 



G(x,x') = (4n)- d/2 A 1/2 (x,x') J dtt~ d/2 exp 



2t 



tiUQ(t\x,x'), (59) 
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where tt(t\x,x') is the transfer function of the operator F (|54|). It is not difficult to show 
that due to the presence of the projector II and the special form of the potential term Q 
(|55|) the asymptotic expansion of the function tr Iiri (t) as t — > begins at the order t N ~ l . 
Comparing the eq. (|5^) with eq. (|T5| ) we find then the function A(t|x, x') in terms of the 
transfer function Q(t\x, x') 

A(t) = (n- i)\r N+1 trun(t). (60) 

By using the expansions of the functions A(t) (|2T| ) and fl(t) ( |T8D and equating the powers 
of t we obtain finally the coefficients c k in terms of the HMDS-coefficients 

c k (x,x) = / fc+jV-1 x tr (n b k+N -!(x, x')), (61) 
\ k J 

where b n are the HMDS-coefficients of the Laplace type operator F (^) with the potential 
term Q (|55|). 

This formula gives the coefficients of the asymptotic expansion of the function y4(t|x, x') 
in terms of well known heat kernel coefficients of the Laplace type operators and deter- 
mines finally the structure of the diagonal singularities of the product of Laplace type 
operators (|28|) . 

From (|29"D we obtain, in particular, the function x(x, x') determining the logarithmic 
singularity 

x ( XjX ') = { - {N ' iy '^m N+lh T (nfy2-i(s,s / )), for even d ^ 
], for odd d 
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